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We present population synthesis modeling of the X-ray background with genetic algorithm - based optimization 
method. In our models the best fit could be achieved for lower values of high-energy exponential cut-off (w 
170 keV) and larger amount of the highly obscured (log Nh ~ 25.5) AGNs. 



1. Introduction 

The cosmic X-ray background (XRB) above 
~ 1 keV is known to be produced by integrated 
emission of discrete sources |26I13I24] . The re- 
sent best reviews in this area are ffl and |14| . 
XRB synthesis models are usually based on the 
so-called unification scheme for AGN 1 , which 
explains the different observational appearances 
by the orientation of accretion disk and molecular 
torus surrounding the nucleus. The intersection 
of the line of sight with the torus determines a 
type 2 AGNs, and the direct observation of the 
nucleus identifies a type 1 AGNs. 

The last works in the area of population syn- 
thesis modeling used assumption about extra 
quantity of absorbed AGNs at high redshifts 
|3()I9I28I1()I2'^ . These models give good approach 
to the exist observations in the energy range of 
1 ^ 100 kev. 

However, there are many unresolved yet prob- 
lems exist, like as a role of the soft X-ray excess in 
AGN spectra at high redshifts, behaviour of the 
AGN luminosity function IT, possible flattening 
of the spectra of AGNs at high redshifts and a 
role of the high luminosity obscured AGNs. Fur- 
thermore, we do not have yet detailed information 
about an exponential cut-off at high energies. 

Here we present population synthesis models 
of the X-ray background in order to investigate 
some additional conditions to obtain the best fit. 



2. Fitting method 

In our work we have tested one of the newest 
effective fitting techniques: a genetic algorithm - 
based optimization 3 . This method implements 
the Darvin's natural selection law for the mathe- 
matical problems. 

The main steps of this technique are: 
1. Constructing a random initial population of 
the sets of the model parameters and evaluating 
the fitness of its members. 2. Constructing a 
new population by breeding selected individuals 
from the old population. 3. Evaluating the fit- 
ness of each member of the new population. 4. 
Replacing the old population by the new popula- 
tion. 5. Test convergence: unless fittest pheno- 
type matches target phenotype within tolerance, 
goto 2. 

The breeding consists of the several steps such 
as crossover (as a result the offspring obtains the 
properties of the both its parents) and mutation 
(which allows to probe the alternative sets of the 
model parameters. This option especially impor- 
tant, when population members become practi- 
cally identical). 

This method already was implemented as one 
of the XSpec fitting methods and now often used 
for different complicated fitting problems. 
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3. Description of the model 

As for any population synthesis model, tlie re- 
sulting spectrum was calculated as a mix of the 
AGN spectra, which are typical for the different 
classes of AGNs, integrated over redshift and lu- 
minosity. To avoid a contribution of the Galactic 
diffuse radiation component, we have modeled the 
energy range above 1.5 keV only. 

Thus, following j4j, for AGN spectra we have 
the next expressions: 

F{Quasars 1) oc * exp( ) (1) 

Ec 

FiSeyferts 1) (x ^-"-^ H<exp(-^) + F,e/K-B)(2) 

F{Type 2) cx F{Type 1) * exp(-crij * Nh) (3) 

where the hard energy index ah = 0.9. For the 
exponential cut-off as a base value we have used 
Ec = 320 keV, but in different models we have 
tried another values. 

The term Frefi [E) represents the Compton re- 
flection component by the accretion disk and has 
been evaluated following ,17j assuming inclina- 
tion angle of 60°. 

For obscured (type 2) AGNs, photoelectric ab- 
sorption cross sections for given hydrogen column 
Nh were calculated following [^Hl- As the base 
distribution of equivalent hydrogen column den- 
sities Nh was taken from 29 (RMS hereafter). 
Following JS| J the local ratio of absorbed and un- 
absorbed AGNs is i? = 4.0 ±0.9. To evaluate this 
ratio with redshift, we have used the next formula 



i.(.)^(^i?(0) + Ai?(oo)*(^l-^^jjx 

X exp(-z/ze) 

which in contrast to the "power law & constant" 
form of |10j has continuous first derivative and 
is convenient to represent the ratio dependence 
from the redshift. Here (i?(0) -f Ai?(oo)) is the 
ratio at "infinity" and distance for R(z) 

to get value i?(0) + ^M2£i. is an exponential 
cut-off scale, the same as in |28|. 



During integration for the X-ray luminosity 
function we have used expression in the form pre- 
sented by 21 . For all our models the Hubble 
constant given by Hf) = 50 km Mpc~^. 

We also computed contribution of the clusters 
of galaxies to the overall background spectrum, 
adopting the temperature distribution taken from 
j2i and cluster luminosity function of 6 . 

Galactic hydrogen column was taken as N^ = 

To fit our models we have used kindly provided 
data from different missions like as HEAO-1 A4 
LED and A2 HED detectors |11I15I2) . ROSAT 
PSPC and ASCA GIS 22.2T]| . 

The main parameters which describe our imple- 
mentation of the genetic algorithm are population 
number (100 members for the model with 3 free 
parameters and 200 members for the case of 4 free 
parameters) , variable mutation rate (initial value 
is 0.005) and high selection pressure regime for 
a breeding of the population members. Further- 
more, we have used so-called " elitism" , when the 
fittest generation member is copied without alter- 
ation into the next generation. A good algorithm 
convergence we achieve after 40-50 iterations. 

4. Results and discussion 

In the present work we realized several kinds of 
models. A simple model A was calculated using 
equation 4 with Ai?(cx)), Zr and Zg as free param- 
eters. The best fit results presented in the Table 
1. The corresponded XRB fit and R{z) curve one 
can see in Fig.l and Fig. 2. 

For the model A we have the same rapid in- 
crease of the R{z) with redshift, as in work of 
|1U| and significantly much more slower exponen- 
tial decay than found by We consider that 
such rapid increase of R{z) can not be real be- 
cause it makes an observer as a particular person. 
The most probable mistake could be using 'as is' 
the distribution of equivalent hydrogen column 
densities of RMS. In spite of very detailed calcu- 
lations of the selection effects, their sample based 
on optical identifications using emission lines and, 
hence, can lose some highly obscured AGNs. In 
reality the distribution could shifts towards high 
values of hydrogen column densities. 
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In order to investigate this possibility in our 
model B we have used as an additional parameter 
the quantity of AGNs with logN^ — 25.5. In the 
Table 1 the last parameter 7125.5 was calculated 
as a fraction of all population of type 2 AGNs 
in units of the RMS fractional population value 
(they have used ri25.5 ~ 0.26). Then new local 
ratio of absorbed and unabsorbed AGNs is calcu- 



lated as 4.0 X (1 



'25.5^ 



'25.5 



) . You can see that 



our new value of 7125.5 is distinctly higher than 
used before. For new model values the XRB spec- 
trum is better fitted below 50 keV but is worse 
above 50 keV. At the same time the R(z) func- 
tion becomes more smooth. The distribution of 
hydrogen column densities for the different cases 
one can see in the Fig. 3. 

The main problem of our fitting curves is high 
bias to the data above 50 keV. In reality if we 
analyse the data up to 400 keV it is obvious that 
we do not have clear exponential cut-off profile. 
Assuming that data below 100 keV represents the 
'true' exponential cut-off we have tried to make 
the model C with Ec ~ 160 keV. One can see 
that in this case we achieve practically ideal ap- 
proach to the observations, although it demands 
to increase the 7125.5 parameter up to 0.93. 

Table 1. The best fit parameters for different 
models. 



model 


AR{oo) 




Ze 


n25.5 


A 


3.34 


0.11 


2.49 




B 


1.21 


3.01 


4.46 


0.47 


C 


2.64 


0.47 


2.96 


0.93 



In spite of approximate character of our mod- 
els we have to conclude the presence of a large 
amount of undetected highly obscured AGNs 
(with log Nh — 25.5). I our future models we 
are going to use more precise simulations of the 
spectra of absorbed AGNs and the latest results 
obtained for the properties of the luminosity 
function. 
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Figure 1. XRB data from different missions and our model fits. 
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Figure 2. The ratio of absorbed and unabsorbed AGNs for different redshifts. 
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Figure 3. The distribution of the hydrogen eohinins of absorbed AGNs for different models and derived 
by [29]. 



